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Abstract 

The two-dimensional cubic nonlinear Schrodinger equation admits a large family of one- 
dimensional bounded traveling- wave solutions. All such solutions may be written in terms of 
an amplitude and a phase. Solutions with piecewise constant phase have been well studied 
previously. Some of these solutions were found to be stable with respect to one-dimensional 
perturbations. No such solutions are stable with respect to two-dimensional perturbations. 
Here we consider stability of the larger class of solutions whose phase is dependent on the spatial 
dimension of the one-dimensional wave form. We study the spectral stability of such nontrivial- 
phase solutions numerically, using Hill's method. We present evidence which suggests that 
all such nontrivial-phase solutions are unstable with respect to both one- and two-dimensional 
perturbations. Instability occurs in all cases: for both the elliptic and hyperbolic nonlinear 
Schrodinger equations, and in the focusing and defocusing case. 

1 Introduction 

The cubic nonlinear Schrodinger (NLS) equation in two spatial dimensions is given by 

^'^pt+Oi^pxx+ Pi'yy + lipl'^ip = 0- (1) 

The NLS equation is said be focusing or attractive in the x-dimension if a > 0. If a < 0, the NLS 
equation is said to be defocusing or repulsive in the x-dimension. Similarly, the sign of /3 leads to 
focusing or defocusing in the y-dimension. The NLS equation is cahed hyperbolic if q/3 < and 
elliptic if af3 > 0. 

Equation (Q) admits a large family of one-dimensional bounded traveling-wave solutions. All 
such solutions, up to Lie group symmetries may be written in the form 

tP{x,t) = (l){x)e'^^^^+'^\ (2) 

where 4>{x) and 9{x) are real-valued functions, and A is a real constant. Solutions of the form ((2)) 
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where c is a real constant. Here k G [0, 1] is the elliptic modulus of the Jacobi elliptic sine function, 
sn{x,k). The function sn (x, A;) is periodic if A; G [0,1), with period given by L = 4i^, where 
K = K{k) is defined by 



and is known as the complete elliptic integral of the first kind. When A; = 0, sn (x, 0) = sin(a;) 
with L = 271. As k approaches 1, sn(a;,A;) approaches tanh(x) and L approaches infinity j^- 
Although (j){x) inherits the periodicity of sn(x,A;), the solution tlj{x,t) is typically not L-periodic 
in the x-dimension, because the periods of e*^ and are typically non-commensurate. 

The solution ip is said to have trivial-phase (TP) if 6{x) is (piecewise) constant and nontrivial- 
phase (NTP) if 6(x) is not constant. Equivalently, the solution ip has TP if c = 0, and has NTP if 
c 7^ 0. For every choice of a and (3, @ specifies a two-parameter family of NLS solutions with the 
free parameters k and B. Without loss of generality, we choose both a and /3 to be ±1. The phase 
contribution 9{x) given in (|3b|) implicitly depends on a and B in both (|Haj) and (|3d|) . In order for 
(j) and 6 to be real-valued functions, we need B G [2A;^,2] if a = 1 or i3 < if q = —1. Figure ^ 
displays the regions of {k, i?)-parameter space that correspond to NTP solutions of the elliptic and 
hyperbolic NLS equations. By varying B so that c — > 0, approaches a (piecewise) constant, and 
the NTP solutions reduce to one of five types of TP solutions: (i) a Stokes' plane wave, (ii) a 
cn-type solution, (Hi) a dn-type solution, (iv) an sn-type solution, (v) a soliton-type solution. The 
limiting solutions correspond to the boundaries of the regions in Fig. ^ Tabled provides the values 
of k and B that cause © to reduce to TP solutions, and also gives the explicit expression for -0 
in the TP limit. Gray solitons, all of which are NTP solutions, result when k = 1 and B < 0. An 
overview of stationary NLS solutions is given in [Hill]. Details of the Jacobi elliptic functions cn, 
dn and sn may be found in 

The stability analysis of TP solutions is well investigated, see for example PUIITTI [TBI ITUl [THI I5] . 
While some TP solutions are stable under one-dimensional perturbations (the bright soliton |15l l9] 
and sn-type solutions jSl), all TP solutions are known to be unstable under two-dimensional pertur- 
bations 0. We know of only the 1981 work of Infeld and Ziemkiewicz 12 for results regarding the 
stability of some NTP solutions of the NLS equation. An additional damping term may lead to sta- 
ble NTP solutions. However, without this term they found that all NTP solutions they considered 
are unstable. Our results agree with their conclusions. The work presented here differs from ^21 in 
that we consider the entire parameter space of NTP solutions, so that all one-dimensional stationary 
solutions of the NLS equation have now been investigated. 

The {k, 5)-parameter space of the NTP solution is two-dimensional, whereas the {k, i?)-parameter 
space of TP solutions is essentially one-dimensional; it forms the boundary of the NTP (k, B)- 
parameter spaces shown in Fig. ^ Thus the TP solutions are only a co-dimension one subset of all 
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(a) a = +1 (b) a = —1 

Figure 1: Admissible parameter space for solutions of the form given in ((2)) for (a) focusing and 
(b) defocusing regimes. The shaded interior region corresponds to NTP solutions. The region in 
the defocusing regime is unbounded below. 



bounded traveling- wave solutions of . The aim of this work is to investigate if has any one- 
dimensional traveling-wave solutions that are stable with respect to either one- or two-dimensional 
perturbations. Since no solutions are known that are stable with respect to two-dimensional per- 
turbations in the TP setting, we focus especially on such perturbations, although one-dimensional 
perturbations are also considered as a special case. 

We investigate the spectral stability of all NTP solutions 0, for all choices of a = ±1 and 
P = ±1. Although all NTP solutions are found to be unstable, our investigations do produce 
important information about the nature of the instabilities of these NTP solutions. 



2 The linear stability problem 

In order to study the linear stability of NTP solutions of the NLS equation, we consider perturbed 
solutions of the form 

^^J^{x, y, t) = (0(x) + eu{x, y, t) + iev{x, y, t) + 0{e^)) e*^(-)+*^*, (5) 

where u{x, y, t) and u (x, y, t) are real-valued functions, e is a small real parameter and (t){x) e«^(^)+*-^* 
is an NTP solution of NLS. Substituting © in linearizing and separating real and imaginary 
parts leads to 

An — 37(/)^n — (iuyy + ac^—ju — 2ac—^4>xV + 2ac—^Vx — auxx = —vt, (6a) 



Xv — j(fPv — Pvyy + ac^—jv + 2ac—7T(l)xU — 2ac—7^Ux — oiVxx = ut- (6b) 
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k value B value 


Solution type 


V' 


a = 1 


k = Be [2/c2,2] 
ke{0,l) B = 2/c2 
k£{0,l) B = 2 
k=l B=2 


Stokes' plane wave 
cn-type 
dn-type 
bright soliton 


VB e*^* 

V2kcn {x,k) e*(2fc'-i)t 
^2 dn {x,k)e'^^-^ ^^ 
\/2 sech(x)e** 


a = —1 


k = B <0 
fee (0,1) B = 
k=l B=0 


Stokes' plane wave 
sn-type 
dark soliton 


^/Be'^^^ 

V2ksn (x, fc)e*(i+^')* 
V2 tanh(x) e*^* 



Table 1: Parameter values {k,B) which reduce NTP solutions to TP solutions. 



Since © does not depend on y or t explicitly, we may assume that u{x,y,t) and v{x,y,t) have 
the forms 

u{x, y, t) = U{x, p, n) e'py+^^ + c.c, (7a) 

v{x, y, t) = V{x, p, n) e'py^^^ + c.c, (7b) 

where p is a real constant, U{x) and V{x) are complex- valued functions, is a complex constant 
and c.c. denotes complex conjugate. Notice that p is the transverse wavenumber of the perturbation 
and O is the exponential growth rate associated with p. If bounded U, V exist such that has 
a positive real part, then the amplitudes of the perturbations grow exponentially in time and the 
unperturbed solution is unstable. 

Upon substitution, © yields the spectral problem 

XU - 3-f(t>^U + f3p^U + ac^-^U - 2ac-^(l)^V + 2ac^V^ - aU^^x = -^V, (8a) 

\V - 7(/>V + + ac^^V + 2ac^ct):,U - 2ac^U^ - aV^^ = nU. (8b) 

If c = 0, then (jSJ reduces to the stability problem for TP solutions. This case is examined in |171 
El n El El El and others. Using the linear system Q, we are now able to investigate the stability 
of all NTP solution numerically. We only consider the stability of NTP solutions in this paper. 
The reader may wish to consult some of the above mentioned references for the stability analysis 
of their limiting special cases. 



3 The numerical investigation of spectral stability: Hill's method 

The main difficulty for the numerical investigation of © is the size of the parameter space in- 
volved. For every choice of the equations parameters a, (5 and solution parameter pairs {k,B), 
the spectrum of needs to be computed for a range of p values in order to determine stability 
or to analyze any instabilities. An efficient numerical method is necessary. Hill's method allows 
for the systematic and efficient exploration of the large phase space encountered here, due to its 
exponential convergence |S1. 
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To apply Hill's method, Fourier expansions are needed for all coefficient functions of (jH)). Using 
the complex Fourier form, we have 



n=— oo n=— oo /q\ 

n.=— oo n=— oo 

where Qn,Rn,Sn and T„ are Fourier coefficients. Note that (j)'^{x) has period L/2 and that is 
never zero except in some TP limit cases. 

The periodicity of the coefficient functions in Q allows us to decompose the eigenfunction 
components U and ^ of (jSJ in a Fourier-Floquet form 

oo oo 

U{x) = e'>''' Y f^ne^"^^/^ and T/(x) = e*/^^ Ke^"^^/^. (10) 

n=— oo n=— oo 

The form of U and V in H1U|) follows from Floquet's theorem and the observation that eigenfunctions 
are bounded, by definition. This decomposition has the benefit of admitting both L-periodic and 
L-anti-periodic eigenfunctions when fj, = 0. Recall that ip is typically only quasiperiodic. Allowing 
fi to be different from gives rise to solutions that are either quasiperiodic or have period greater 
than 2L. 

Substitution of © and (jTU]) into ((SJ and equating Fourier coefficients allows us to write equations 
for Un and Vn as a coupled bi-infinite system of difference equations given by 

(/ imT \ 2\ °° °° 

A + — a lifi-\ — j \ Un + 3j Y^ Q n~m Um. — Oi(? Y^ S n-m U^a 

^ ^ / m=— oo m=— oo 

oo ^ . \ oo 

+ 2qC Y^ Tn-m Vm. — 2aC iifJ.-\ —j Y^ R n-m Vm = ^Vn , (Ha) 

m=— oo ^ ^ m=— oo 



A + Pp"^ - aiifi + -J- j j K - 7 ^ Qn^Vm + ac^ ^ 5 n-m Ki 

s, ^ ^ / m=— oo m=— oo 

oo y . \ oo 

+ 2ac Y^ Tn-mUm — 2ac iifi -\ — —j R n~m Um. = ^Un, (Hb) 

m=— oo ^ m=—oo 



for all integers n. Here /i € [7^;-^) and Q n-m = if ^^^^y^ Z, with i?(.),5(.) and T(.') defined 
similarly. The system of equations is equivalent to the original system (jH)) . 

Remarks 

• In practice, a pre-multiplication of the linear system by (j)^ allows for the exact Fourier series 
expansion of (/)^, cf)^ and (j)^ to be used. This follows from the differential equations for sn (x, k) 
and Jacobi's series expansion of sn ^(x, k) This pre-multiplication transforms the original 
eigenvalue problem into a generalized eigenvalue problem jllj . 
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Parameter 


Description 


Value 


k 


Elliptic modulus 


linspace (0 , 1 , 65) 


B 


bniit 


For a = —1 : — logspace(— 8, 0, 65) 

For a = 1 : {2k'^ + logspace(-8, 0, 65)) n {2k'^ , 2) 


N 


Fourier cutoff 


For a = -1 : 15 + ceil(5k^) 
For a = 1 : 10 + ceil(25ki°) 


P 


Perturbation wavenumber 


linspace (0 , 4 , 65) 


P 


Floquet parameter 


linspace(-f,f,2l) 



Table 2: Parameter values and ranges used in numerical experiments. Only perturbations of NTP 
solutions are considered. 

• Note that Hill's methods enables one to compute the spectrum of a linear operator with 
periodic coefficients. Despite the fact that the solution ^ is typically quasiperiodic, the 
coefficient functions of the linear stability problem (jS)) are always L/2-periodic. 

4 Numerical experiments 

By choosing a finite number of Fourier modes, the exact bi-infinite system ()11() is truncated. We 
explicitly construct and compute approximations to the spectral elements (i.e. eigenvalues or 
elements of the continous spectrum) of by finding the eigenvalues of the truncation of (|1H) . We 
consider all four cases individually: (I) focusing in both x and y {a = (3 = 1), (II) focusing in x and 
defocusing m. y {a = —j3 = 1), (III) defocusing in x and focusing in y (—a = /3 = 1) and finally, 
(IV) defocusing in both x and y (—a = —j3 = 1). 

In each case, a large number of parameter values in the two-dimensional parameter space shown 
in Fig. n was explored numerically. The (/c, i?)-parameter values considered correspond to NTP 
solutions, and do not include TP solutions or gray solitons. Approximately 5.2 million generalized 
eigenvalue problems were considered, the size of each determined by the cutoff mode N of the 
underlying Fourier series. A truncation to N positive Fourier modes reduces the exact bi-infinite 
system (jllj) to an (4A^ -|- 2)-dimensional approximate problem. For several choices of k and B, a 
value of = A^(A:, B) was chosen to ensure that the resulting eigenvalues had converged to within 
a measured tolerance. A simple polynomial was used to fit this data. This information, and details 
related to other parameter ranges used in the experiments, are included in Table [2 In the table, 
k is the elliptic modulus and B is the offset parameter (as in ©), (4 A -|- 2) is the dimension used 
to approximate (jH)), p is the wavenumber of the perturbation in the y-dimension, and /i is the 
Floquet exponent. Lastly, linspace(a, b,m) is a linearly spaced vector from a to 6 of length m, 
logspace(a, b,in) is a logarithmically spaced vector from 10" to 10^ of length m and ceil(x) is the 
smallest integer greater than or equal to x. 

5 Observations from the numerical investigation 

First and foremost, it should be stated that none of the NTP solutions considered here were found to 
be spectrally stable under one-dimensional (p = 0) or two-dimensional {p ^ 0) perturbations. This 
establishes, at least numerically, that all one-dimensional traveling-wave solutions of NLS of the 
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form are spectrally unstable with respect to either one- or two-dimensional perturbations. At 
this point, it remains to investigate the nature of the instabilities and their corresponding growth 
rates, so as to better understand the dynamics of this important class of solutions of the NLS 
equation. 

Using Hill's method we numerically considered the instabilities due to perturbations with 
wavenumber p € [0,4]. Note that p = corresponds to one-dimensional perturbations. For each 
NLS equation (i.e. for each choice of a,/3 = ±1), and solution (i.e. each parameter pair {k,B)), 
and for each perturbation of wavenumber p, an equally-spaced sequence of Floquet parameters p, 
was chosen from the interval [— ;^]- The generalized eigenvalues and eigenvectors were computed 
for the matrix that results from a truncation of The generalized eigenvalues are approxima- 

tions of spectral elements of 0, and an approximation of the corresponding eigenfunctions may 
be reconstructed from the eigenvectors. 

Since a single eigenvalue with positive real part leads to instability of the system, the approxi- 
mate eigenvalue with largest real part over all choices of p was recorded for each (k, B, p) triplet. 
That is, we compute 

f^growth(A;,S,p) = max Re[0(A;,B,p,/i)] , (12) 

^l&[—^r/K,■K/K] 

which we call the (most unstable) growth rate, which we call the (most unstable) growth rate. 
The value Ingrowth represents the largest exponential growth rate a given NTP solution with pa- 
rameters (/c, B) will experience when perturbed with transverse wavenumber p. It also allows us to 
determine the perturbation to which the NTP solution is spectrally the most unstable. We reduce 
the dimension further by computing the largest such growth rate over all sampled perturbation 
wavenumbers p. This quantity, 

^ma.yi{k,B) = max Ogi.owth(/c, -B, /o) , (13) 

pe[o,4] 

the maximal growth rate over all p, is plotted in the first columns of Figs. |21 and |21 The value 
^max represents the maximal exponential growth rate that a solution with parameters {k, B) can 
undergo in the range examined. We also recorded the minimum growth rate over all />, 

Q.rain{k,B)= min Ogrowth(A;, -B, p), (14) 
pe[o,4] 

to verify that all solutions are unstable with respect to every sampled perturbation. 

Every point plotted in Figs. [21 and corresponds to an NLS solution for which we considered 
the linear stability analysis, and the boundaries in the figures are the boundaries of the regions 
represented in Fig. ^ corresponding to limiting TP solutions. Fig. [21 corresponds to the x- focusing 
(a = 1) parameter range {k,B) = (0,1) x (2A;^,2) in the q = 1 case of Fig. The one-to-one 
transform Tf{B) = (B — 2/c^)/(2 — 2/c^) is used to normalize the range of B. This maps the 
interval [2A;^,2] to [0,1]. Figure [31 corresponds to the x-defocusing (a = —1) parameter range of 
{k,B) = (0,1) X (-1,0) shown in Fig. [Ha. The transform Td{B) = -B is used in Fig. [3 A logio 
scale is used in the vertical dimension of Figs. [21 and [21 This causes the panels of Fig. [21 to become 
increasing sparse in their lower right corners. The right-hand panels of Figs. [21 and [21 indicate 
the wavenumber p that leads to maximal growth shown in the left-hand panels. Recall that our 
computations were limited to p E [0, 4]. 
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5.1 Case I: Elliptic setting with a = P = 1 

Panels la and lb of Fig. [2 summarize some properties of the computed instabilities in the case 
of focusing in both the x- and y-dimensions. The lower boundary of the plot corresponds to 
B = 2k'^ + (10~®), and is therefore only slightly away (in the parameter space of B) from a cn-type 
solution. The upper boundary is close to dn-type solutions, with B = 1.99. The left boundary of 
the plots, where k = 0.01, represents a region in parameter space near to Stokes' wave solutions. 
The entire right-hand boundary, where k = 0.99, is near to the bright soliton limit case which 
occurs at {k,B) = (1,2). 

A distinct ridge of large instability is noticeable in the plot of J^max in panel la of Fig. ^ The 
ridge appears to begin near the zero solution at {k,B) = (0,0), and remains close to the cn limit 
boundary (within approximately .02 units, remembering the log^^g scaling) as k increases. Moving 
away from the cn boundary results in the rapid increase of ^^max- Movement away from the dn 
boundary results in a slower increase in the value of J^max) as does moving away from the Stokes' 
wave boundary for B greater than approximately 0.001. The maximum value of f^max over the 
sampled {k,B) space, given by iimax = 5.666, is reached at {k,B) = (0.99,1.98). This growth 
rate should be compared to the maximal growth rate of the corresponding TP case [5| which is 
i?max = 1- The minimum (iJ^in = 0.015693) occurs for {k,B) = (0.01,0.01) for p = 4. Note that 
this is on the boundary of the computational domain and that at (k, B) = (0, 0), corresponding to 
the spectrally stable = solution, i?min = 0. 

In panel lb, the wavelength corresponding to the maximal growth of la is given. In this case, 
the maximum instability occurs for the shortest wavelength samples, p = A. This indicates that 
there is a strong short-wavelength instability. 

5.2 Case II: Hyperbolic setting with a = —[3 = 1 

Panels Ila and lib of Fig. El summarize some properties of the computed instabilities in the case 
of focusing in the x-dimension and defocusing in the y-dimension. The lower boundary of the plot 
corresponds to -B = 2fc^ -|- (10~^), and is therefore only slightly away (in the parameter space of B) 
from a cn-type solution. The upper boundary is close to dn-type solutions, with B = 1.99. The 
left boundary of the plots, where k = 0.01, represents a region in parameter space near to Stokes' 
wave solutions. The entire right-hand boundary, where k = 0.99, is near to the bright soliton limit 
case which occurs at {k,B) = (1,2). 

As in Case I, a ridge of large growth rate is noticeable in the growth plot shown in panel 
Ila. The ridge appears to begin near the zero solution at {k,B) = (0,0), and remains close to 
the cn-type limit boundary (within approximately .02 units, remembering the log^g scaling) as k 
increases. This ridge has a local minimum near k = 0.7 and increases to a global (over all admissible 
{k, i?)-parameter space) maximum at A; = 0.96. As in the setting above, moving away from the 
cn-type boundary results in a rapid increase of rJmax- Moving away from the dn boundary results 
in a slower increase in the value of i^max- For B > 0.001, moving away from the boundary result 
in a similar slow increase in Jlmax- For k > 0.96, it appears that the limiting value of f^max is 
consistent with the bright soliton results of |151 19]. The maximum (-Rmax = 6.1141) and minimum 
(Rmin = 0.012535) growth rates span a slightly larger range than the similar values in Fig. la. 
These occur at {k,B) = (0.01,0.01) and {k,B) = (0.96,1.98), respectively. The maximal growth 
rate should be compared to the maximal growth rate of the corresponding TP case 5^ which is 

-^max — !■ 
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la. r^max with (a,/3) = (1,+1) lb. Corresponding = (1,+1) 



R . =0.015691 

mm 




Wave number 1 1 




p = 4 



01 02 03 04 05 06 07 08 09 



01 02 03 04 05 06 07 08 09 



Ila. Omax with (a,/3) = (1,-1) lib. Corresponding p, {a, (3) = (1,-1) 




Wave number 1 -1 




01 02 03 04 05 06 07 08 09 







01 02 03 04 05 06 07 08 09 



Fi gure 2: Focusing in the x-dimension. The first column contains surface plots of r^max vs. {k, B) 
and the second column contains the surface plots of the maximizing wavenumber p vs. {k,B). 

R-max — m&X/j R ^^max and Rmin 

= miuyt^s ^^min- White space corresponds to (A;, i?)-parameter space 

that was not sampled. 



In panel lib, the wavelength corresponding to the maximal growth i?max of Ila is given. In this 
case, the maximum instability of Rmax occurs for p = 3.375. The surface shown in lib appears to 
be more smooth than the surface of lb. 



5.3 Case III: Hyperbolic setting with —a = P = 1 

Panels Ilia and Illb of Fig. IHl summarize some properties of the computed instabilities in the 
case of defocusing in the x-dimension and focusing in the y-dimension. The lower limit of the 
plot corresponds to B = — (10~^), and so is just slightly away from the sn-type solution. The left 
boundary of the plots, where k = 0.01, represents a region in parameter space near Stokes' wave 
solutions, while k = 0.99 on the right boundary is near to the gray soliton limit. 

A distinct ridge of large instability is noticeable in the growth plot displayed in panel Ilia. The 
ridge appears to begin near the zero solution limit at {k,B) = (0,0), and remains close to the sn 
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limit boundary (within approximately .02 units, remembering the log^Q scaling) as k increases. It 
quickly reaches the global maximum (ovgi all admissible (A^, -B)-parameter space) of -f^max — 

7.6375 

at k = 0.02 and B = —0.0001. The ridge then decreases in amplitude as k increases towards 1. 
Moving away from the sn-type boundary results in a rapid increase of i^max- Moving away from the 
dn-type boundary results in a slower increase in the value of i^max- Similarly, the increase is slower 
when moving away from the Stokes' wave limit for B < —0.001. The maximum exponential growth 
rate, .Rmax = 7.6375, occurs for {k,B) = (0.02,-0.00001). This growth rate should be compared 
to the maximal growth rate of the corresponding TP case [Jj which is -Rmax = 1- The minimum 
exponential growth, Rmin = 0.015578 is found at {k,B) = (0.01,-0.9). Both the maximum and 
minimum are located near the Stokes' wave boundary. By restricting p = and allowing B to 
approach zero, f^max — *■ 0, and the one-dimensional stability result of the sn-type TP solution of [2] 
is recovered. 

The plot Illb indicates short -wave perturbations lead, to large values of ^^max* 

The largest 

growth occurs for a perturbation with wavenumber of p = 3.625. 
5.4 Case IV: Elliptic setting with —a = —(3 = 1 

Panels IVa and IVb of Fig. 01 summarize some properties of the computed instabilities in the case of 
defocusing in both the x- and y-dimensions. The lower limit of the plot corresponds to B = — 10~^, 
and so is just slightly away from the sn-type solution. The left boundary of the plots, where 
k = 0.01, represents a region in parameter space near to Stokes' wave solutions, while k = 0.99 on 
the right boundary is near to the gray soliton limit. 

A distinct ridge of large instability is noticeable in the growth plot shown in panel IVa. The 
ridge appears to begin near the trivial limit k = and B = 0, and remains close to the sn limit 
boundary (within approximately .02 units, remembering the log^g scaling) as k increases, to reach 
a global maximum at A; = 0.02 and B = —0.0001. The ridge then appears to decrease in amplitude 
as k increases towards 1. As in Case III, moving away from the sn-type boundary results in a 
rapid increase of r^max- Moving away from the dn-type boundary results in a slower increase in 
the value of J^max- The same is true when moving away from the Stokes' boundary, when B is 
less than approximately -0.001. The maximum exponential growth, rate, Rmax — 

7.6456, and the 

minimum, Rmin = 0.0001556, span a slightly larger range of values than do the values of Omax 
in panel Ilia. The maximum and minimum values are obtained at {k,B) = (0.01,-0.00009) and 
{k,B) = (0.01, —1), respectively. Both are located near the Stokes' wave boundary. The maximal 
growth rate should be compared to the maximal growth rate of the corresponding TP case [7] 
which is -Rmax ~ 0.26. As in Case III, restricting p = and allowing B to approach zero results in 
^max 0, and the one-dimensional stability result of ,5 for the sn-type TP solution is recovered. 

In panel IVb, wavenumbers corresponding to f^max of IVa are given. It appears that a majority 
of the large values of Omax are attributable to small-p (long-wave) perturbations. In fact, the 
largest growth occurs for p = 0, the one-dimensional perturbation. This should be contrasted with 
the three previous cases, where short-wavelength two-dimensional perturbations with wavenumber 
p > 3 were associated with the largest Omax values. 
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Ilia. Qmax with {a,p) = (-1,+1) 

ab -1 1 

li 




Illb. Corresponding p, {a, (3) = (—1,+!) 



Wave number -1 1 




01 02 03 04 05 06 07 08 09 



01 02 03 04 05 06 07 08 09 



IVa. Omax with (q,/9) = (—1,-1) IVb. Corresponding p, (a,/3) = (—1,-1) 




I 



Wave number -1 -1 




01 02 03 04 05 06 07 08 09 



01 02 03 04 05 06 07 08 09 



Fi gure 3; Defocusing in the x-dimension. The first column contains surface plots of ^^max ^s. 
{k,B) and the second column contains surface plots of the maximizing wavenumber p vs. {k,B). 

R-max = maXfc^B J^max and Rmin = miUfe^B ^min- 



11 



6 Summary 



In this paper, we considered the spectral instabihty of one-dimensional traveling-wave nontrivial- 
phase (NTP) solutions of the cubic nonlinear Schrodinger equation. Such solutions are expressed 
in terms of Jacobi elliptic functions. An exact spectral form of the linearized operator is truncated 
and used to construct an associated generalized eigenvalue problem. The positive real parts of the 
resulting eigenvalues were used to determine that there are no stable NTP solutions. 

Numerical results indicate a well-defined ridge of large growth rate located in the {k,B)- 
paramctcr region associated with nontrivial-phasc solutions. This implies that the most unstable 
NTP solutions are more unstable than any TP solution, in the sense that they exhibit larger ex- 
ponential growth rates. Further, for all cases the exponential growth rate fimax increases when 
moving away from the limiting TP solutions. This divergence is gradual in some cases, but very 
sharp in other cases, as discussed above. 

In summary, numerical evidence suggests that all bounded, nontrivial-phase one-dimensional 
traveling-wave solutions of the cubic NLS equation are unstable with respect both one-dimensional 
and two-dimensional perturbations. 
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